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Abstract 

We investigate the electronic dispersion of high-T c superconductor on the basis of the two- 
dimensional three-band Hubbard model with the electron-phonon interaction together with the 
strong electron-electron interaction. In our model, it is shown across the hole-doped region of 
high-T c superconductor that the electron-phonon interaction makes a dispersion kink, observed 
along the nodal direction, and that the small isotope effect appears on the electronic dispersion. 
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INTRODUCTION 

For the past two decades, extensive studies of high-T c cuprates have spotlighted many 
curious phenomena. It has been argued that most phenomena are attributable to the strong 
correlations among electrons, which play significant roles in these materials. However, since 
the discovery of sudden changes in the electron dispersion or "kinks" shown by the angle- 
resolved photoemission spectroscopy (ARPES),-^ effects of electron-boson interactions on 
electronic self-energies have been recognized. While the kinks are now indisputable in 
cuprates,- their origin as arising from electronic coupling to phonons 4 or magnetic exci- 
tations^^ remains unclear. 

Recent scanning tunneling microscope study showed that the statistical distribution of 
energy of bosonic modes (Q) has meaningful difference between 16 O and 18 O materials. 
Thus, it should be hard to exclude the possibility that electron-phonon interaction (EPI) 
significantly affects the electronic states in cuprates. 

In this study, we investigate the analysis upon the EPI together with the electron-electron 
interaction (EEI) on the basis of the two-dimensional (2D) three-band Hubbard-Holstein 
(HH) model. With the use of our three-band HH model, we can reproduce the situation of 
the real high-T c materials, in which EPI mainly works on p electrons at O sites. 

FORMULATION 

Our model Hamiltonian H is composed of d electrons at each Cu site, p electrons at O 
site, and lattice vibrations of O atoms. We consider only the on-site Coulomb repulsion U 
between d electrons at each Cu site as our EEI. Let us define that \x and N represent the 
chemical potential and the number of all electrons, respectively. Then, H — fiNo is divided 
into the non-interacting part, Hq, the electron-electron interacting part H e \- e \, the phonon 
part ifph, and the electron-phonon interacting part i^ei-ph as 

H — fiN = H + H e i_ d + H ph + H e i_ ph , 

No = E(4A* + Vtvl a + vtM- (1) 

kcr 

Here, (iko-t^Lr) an d Pkcf (P^a ) are the annihilation (creation) operator for d and p x ^ elec- 
trons of momentum k and spin a, respectively. The non-interacting part Hq is represented 
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We take the lattice constant of the square lattice formed from Cu sites as the unit of length. 
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is the transfer energy 



between a d orbital and a neighboring p x ( y ' orbital and t pp is that between a p x orbital and 
a p y orbital. A^ p is the difference of energy levels of d and p orbitals. In this study, we take 
tdp as the unit of energy. The residual parts are described as follows: 
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where U is the on-site Coulomb repulsion between d orbitals, N is the number of k-space 



lattice points in the first Brillouin zone (FBZ), and g a\, [y = {x, y}) is the electron-phonon 
matrix element, respectively. We consider that the half-breathing phonon mode, 9 in which 
oxygen ions are vibrating along the x or y directions, is crucial for our problem. Thus, 
ignoring the other phonon modes, we have the electron-phonon interacting part as Eq. (J5j). 

Then, we introduced the unperturbed and perturbed Green's functions, which are to be 
described in 3 x 3 matrix form. The unperturbed Green's function G^(k, z) is derived from 
Eq. (|2D as 



G(°>(k, 



\zl - H ( 



i-i 



(6) 



where I is a 3 x 3 unit matrix. Using the abbreviation of Fermion Matsubara frequencies, 
e n = nT(2n + 1) with integer n and temperature T, the perturbed Green function G(k, z) 
is determined by the Dyson equation: 



G(k, ie n ) - 1 = G(°) (k, ie n ) - 1 - S(k, ie n ) 



(7) 
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where S(k, ie n ) is the self-energy expected to be a diagonal matrix. In order to estimate the 
d electron self-energy E^k, ie n ) = £n(k, ie n ) in Eq. ([7]), we adopt the fluctuation exchange 
approximation^ as follows: 

T 

S d (k, ie n ) = — ^ G d (k - q, ie n - iu m ) 14i-d(q, i^m), (8) 



N 

- 1 " qm 



where G d (k,ie n ) = G n (k,ie r 



XT- / • \ 3 f/ 2 x(q,i^m) . 1 £/ 2 x(q,i^ 
l4i_ e i(q, iw m ) = — 7 — : r + - 



2 1 - C/x(q, i^m) 2 1 + C/x(q,io; TO ) 
-f/\(q> m ), (9) 

where io> m = 2m ttT with integer m are Boson Matsubara frequencies, and 

T 

X(q, iu) m ) = -jz G rf(q + k ' ' luJ ra + ie„) G d (k, ie n ). (10) 

k n 

In order to estimate the p^-electron self-energy E x (y) (k, ie n ) = ^22(33)(k, ie n ) in Eq. (J7|), 
we exploit the Brillouin-Wigner perturbation theory. We adopt the self-consistent one-loop 
approximation as follows: 

T 

Ex(y)(k, ie n ) = — G x{y) (k - q, ie n - iu m ) Ki- P h(q, ^m), (H) 



qm 



where GO^k, ie n ) = (^22(33) (k, ie n ) and V^ip h (q, iw m ) is the EPI on p x ^ electron. Our EPI 



c (y)lK, Lt n ) = U- 22 (33)l^,ifcnJ ailU K el _ ph 

is determined as follows 
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where A = g 2 /(2u^) and a^*" = sin-^^-. uj^ is the specific phonon energy for the half- 
breathing mode. As above, we ignore the effects in which EEI and EPI are coupled. Thus, 
as shown by Eqs. ([8"])- (jl2p . in our formulation, EEI and EPI are completely decoupled. Of 
course, this assumption is inadequate to analyze the case in which the characteristic energy 
due to EEI is comparable to the phonon energy. However, as will be seen later, we actually 
treat the cases with rather high characteristic energy due to EEI. Hence, decoupling EEI 
with EPI should be justified in our analysis. The electron-phonon coupling constant A does 
not depend on Mo since 00^ oc Mq 1//2 and g oc (MotUh) -1 ^ 2 , where Mq is the mass of an 
oxygen ion. In our model, thus, the isotope effect is reflected on the phonon energy in the 
electronic self-energies only, but not any changes in the strength. 
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16q 


0.5575 


0.4837 


0.0412 


LD 


18q 


0.5574 


0.4838 


0.0412 


LD 


16q 


0.5932 


0.5039 


0.0971 


UD 


18q 


0.5931 


0.5040 


0.0972 


UD 





0.6281 


0.5258 


0.1540 


OP 


18Q 


0.6281 


0.5259 


0.1540 


OP 


16Q 


0.6683 


0.5525 


0.2209 


OD 


18Q 


0.6683 


0.5526 


0.2209 


OD 


TABLE I: Number 
RESULTS AND DISCUSSION 


of doped holes. 


S = n h A + n£ - 1. 





We need to solve Eqs. (ITl)- (llll) in a fully self-consistent manner. During numerical cal- 
culations, we divide the FBZ into 128 x 128 meshes. We prepare 2 12 = 4096 Matsubara 
frequencies for temperature T ~ 87K. As shown later, at this temperature, our calcula- 
tion can reproduce the important behavior of electrons in normal state. Moreover, to our 
knowledge, the situation will not be changed if we change T to some extent. 

td p ~ 1.0 eV and t pp ~ 0.55 eV, which are all common for our calculations. These values 
are chosen so that we can reproduce the typical Fermi surface of Bi 2 Sr 2 CaCu 2 08+ ( 5 observed 
by ARPES.^2 A dp ~ 1.4 eV, U ~ 3.0 eV, and A = 0.8 unless stated. The phonon energy is 
set as ~ 65(61) meV for 16 0( 18 0) material. 

We show the numbers of doped holes for our fully self-consistent solutions in Table [H The 
numbers of doped holes both for 16 and 18 materials are exactly the same to three places 
of decimals and they correspond to four different hole-doped samples, lightly doped (LD), 
underdoped (UD), optimally doped (OP), and overdoped (OD), respectively. In Fig. [TJ we 
show the color map of the one-particle spectrum at Fermi level A(k, 0), where 

A(k, e ) = -ilm{TrG(k,i6 n )} i „, (13) 

in order to indicate the Fermi surfaces for 16 materials. In Eq. (TT3]) . the Pade approximation 
is exploited for analytic continuation. Furthermore, we calculate the electronic dispersions of 
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FIG. 1: Fermi surfaces for 16 materials indicated in Tabled! 

the antibonding band E^(k) along the nodal direction indicated as the cut in Fig.[T]for all our 
doping cases. Ed(k) is determined as the k point on which A(h, e) has the maximum value at 
each energy level. Due to this method, the curves of £/rf(k) look like a series of line segments, 
as shown in Fig. [21 We compare every i^(k) of our solution with the one obtained by another 
fully self-consistent manner, in which the same calculation is performed, except for the EPI. 
We show our results on these dispersions in Fig. [2j where we can easily recognize that the 
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-0.08 -0.06 -0.04 -0.02 0.00 
k-k F in (jt,7t)/a units 

FIG. 2: Dispersion kinks along the nodal direction. Momentum is measured from each Fermi 
surface. Arrows indicate the momenta at which kinks occur. 
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FIG. 3: Doping dependence of kink energies. Bars indicate the discretization error during analytic 
continuation. 




-0.08 -0.06 -0.04 -0.02 0.00 
k-k F in (7t,7t)/a units 

FIG. 4: Dispersion kinks along the nodal direction for A = 1.2. The number of doped holes for 
16 0( 18 0) material is 5 = 0.1309(0.1310). The dashed lines show the ones for OP in Table QJ 

dispersion kinks along the nodal direction appear only when EPI affects the p x ^ electrons. 
The kink energies were slightly shifted by 16 — > 18 substitution. In Fig. [3], we detail how 
these kink energies shift depending on hole doping. These theoretically evaluated isotope 
shifts are at most 2.5 meV, which are much smaller than the ones measured by another 
group's ARPES experiment.— 1 ^ Furthermore, these isotope shifts are almost independent 
of hole doping while another group insists that they are critically affected.— Considering the 
energy and momentum resolutions in their experiment, it may be hard to detect the subtle 
isotope shifts and their dependence on hole doping shown in our model. 

Let us now look at A and dependences in the dispersion kinks along the nodal direction 
in detail. Figure H] shows the energy dispersions for A = 0.8 (dotted lines) and A = 1.2 (solid 
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FIG. 5: Dispersion kinks along the nodal direction for A^„ = 1.8 eV. The number of doped holes 
for 16 0( 18 0) material is S = 0.1621(0.1622). The dashed lines show the ones for OP in Table E 
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FIG. 6: E p (k, u) on k = (%p — j^j (jr, tt) (n = 0,1,..., 5), where % = ^ for (a), fcp = ]§g for 



(b), fcp = ^ for (c), and &f = t§5 for (d), respectively. 
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FIG. 7: S<i(k, u) on the same k points as in Fig. |U 
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FIG. 8: Dispersion kinks along the nodal direction for U = 2.5, 3.0, and 3.5 eV. They all correspond 
to OP. Momentum is measured from the Fermi surface for U = 3.0 eV. 
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lines). There is no clear difference in the isotope effect between A = 1.2 (2.5 meV) and 
A = 0.8 (2.0 meV), though the dispersion kinks for A = 1.2 are distinctly shifted to the 
high binding energy compared with those for A = 0.8. On the other hand, the isotope shift 
for A^p = 1.8 eV (lmeV) is slightly shrank compared to that for A^ p = 1.4 eV (2.5 meV), 
though the dispersions depend on the A dp considerably, as shown in Fig. Hence, our 
model calculation shows that the isotope shifts are not sensitive to A and Ad p . Considering 
that Adp is closely related with EEI in our three band HH model as discussed later, we can 
be fairly certain that the isotope shifts are determined by the relative strength between EPI 
and EEI. However, these changes of the isotope shifts are minute, thus, our discussions so 
far are valid regardless of A and A dp . 

To clarify the EPI effect on the p x ^ electrons described above, we investigate the p 
electron self-energy E p (k, u) = E-^k, uS) + E 2/ (k, uS) along the nodal direction. In Fig. [61 we 
show E p (k, lo) on every six k points along the nodal direction, located inside Fermi surfaces. 
The energy where S p (k, uj) is maximal corresponds to the one of the dispersion kink and 
shifts upward by 16 O — > 18 O substitution, as shown in Fig. [2j The energy dependence of 
S p (k, uj) is definitely due to the EPI introduced with the use of Eqs. (JTj), ( fill , and ( fT2l) . 
Thus, we can conclude that in our solutions for all doping levels from the UD to the OD 
region, the dispersion kinks along the nodal direction are created only when the EPI is 
included. 

Hereafter, we will discuss why the magnetic ingredients hardly bring the dispersion kinks 
along the nodal direction. Even when EPI does not exist, the electrons in our results are 
exposed to the strong AF fluctuation originating from the electronic correlation among 
d electrons. In Fig. [71 we show E^(k, on the same six k points as in Fig. [6j It is 
shown that there is no clear difference in Ed(k, uj) between 16 and 18 materials. The 
energy dependence of E d (k, uj) is definitely due to the electronic correlation introduced with 
the use of Eqs. (l8l)-(fT0l. We easily recognize that Ed(k, uj) uniformly increases with the 
binding energy and has no maximal value up to 200 meV even for LD case, in which the 
strong AF fluctuation is expected to be grown. Thus, along the nodal direction, the strong 
AF fluctuation could cause the renormalization of the Fermi velocity, however, it hardly 
promotes any anomalous behavior such as kink structure. 

Finally, we will discuss how our EEI affects the dispersion for 16 O materials. When the 
on-site Coulomb repulsion U is changed, the Fermi velocity is renormalized differently, but 
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this would not change the kink energy so much since the kink energy is determined by the 
EPI alone. In Fig. [HI we lay out our results for three different Us and they all correspond 
to OP. Their total doped holes 5 are slightly different as: 5 = 0.158,0.154, and 0.149 for 
U = 2.5,3.0, and 3.5 eV, respectively. Hence, the Fermi momentum for U = 2.5 eV moves 
inside (or the Fermi surface shrinks) and the one for U = 3.5 eV moves outside (or the 
Fermi surface enlarges), compared to the one for U = 3.0 eV. These changes of the Fermi 
momenta are small; however, the dispersions at higher energy are quite affected, reflecting 
the binding energy dependence of S^(k, to), as shown in Fig. [7J Therefore, the dispersion 
at higher energy could be changed a lot by the EEI even if the Fermi velocities are almost 
independent of them. 



CONCLUSIONS 



By the analysis of our model, we can show that the dispersion kink along the nodal direc- 
tion occurs due to EPI. The isotope effect upon the electronic dispersion is shown near the 
kink of energy dispersion, not in the high binding energy portions. — ^ Our evaluation of 
the subtle isotope shifts has been backed by the report of the recent ARPES experiments^ 1 ^ 
which show the lack of the unusual isotope effect in the high energy portion.— Fortu- 
nately for us, our scenario was possibly realized in further ARPES experiment.— In addition 
to that, we have investigated how EEI effects on the nodal dispersion. It can hardly af- 
fect the kink and the nodal Fermi velocity; however, it can change the dispersion at higher 
energy. Hence, EPI and EEI play different roles on the nodal energy dispersion, respectively. 

Of course, our treatment of EEI is just suited for weak coupling regime, and all of our 
parameter sets employed might be far from the ones for strong coupling regime. If we 
investigate strong coupling regime with the use of another approach, the low energy structure 
corresponding to the superexchange J could appear in the dispersion. However, our results 
presented here suggest that the structure should appear as a broad peak at higher energy 
due to the frequency dependence of the strong AF fluctuation, which will grow into J. 
As we all know, other works have already derive qualitatively similar conclusions on the 
basis of other models. Some groups adopt t-J models i£i22i2ii22 anc l other groups do one- 
band HH models.— Mi2k Furthermore, other groups have succeeded in explaining the ARPES 
results.— i2L22i22a2Jli2i However, in our 2D three-band HH model, both the electron-electron 
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interaction among the d electrons and the EPI on p electrons are considered according to 
high-T c materials. We believe that it is important that quantitatively consistent results with 
the ARPES experiments can be reproduced from such a model. The advantage will be when 
our analysis extends to the superconducting state, in which the p electrons play important 
roles as well as d electrons. 
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